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NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 



TECHNICAL MEMORANDUM NO. 8 61 



COMPRESSION STRUTS 



WITH NONPROGRESSIVE!! VARIABLE MOMENT OF INERTIA* 



By 3. Radomski 



SUMMARY 



The buckling failure conditions for a "bar with non- 
pr ogressively variable moment of inertia J n , although 
constant over length l n , are developed. 

For two cases: 1) bar consisting of two lengths l x 
and l a with J i and J 2 ; 2) bar consisting of three 
lengths l lt l s and again l x , with J lt J s and again 
J i (symmetrical with respect to center), graphs are 
plotted for different ratios Jj. : J g over l x : l 2 , show- 
ing a mean moment of inertia J m , with the aid of which 
the buckling strength P^ of the bar with sudden varia- 
tions of the moment of inertia can be represented in the 
Euler form. 



I. INTRODUCTION 



only valid in the 



(fig. 1) 

P k is buckling load at failure, f m in» minimum section. 

In addition, it is assumed that the moments of inertia at 
the points of sudden variations are effective in their 
full magnitude immediately before and after the jump, 
whereas in reality there is a brief compensating zone for 
welded and a longer zone for riveted bars. 



The following investigation is 
Euler range, that is, only when 



For the rest, the assumptions are that the bar was 

♦"Knickstabe mit sprungweise veranderlichem Tragheits- 
moment . " Luf tf ahrtf orschung, vol. 14, no. 9, September 
20, 1937, pp. 438-443. 
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previously exactly straight' and perfectly pin-jointed, 
that the effect of the cross force and the change in 
buckling length due to axial load and floxuro is negli- 
gibly small and that the differential equation of the 
elastic line E J y" = - M is applicable. 



II. GENERAL DERIVATION 



Consider an originally straight bar, consisting of 
n joined parts stiff in bending with n different - but 
along the lengths l n constant - moments of inertia. The 
axial load P has an eccentricity e of finite magni- 
tude, small in relation to the cross-sectional dimensions, 
that has no effect on the result of the investigation and 
is, for the sake of simplicity, constant (fig. 2). 

Jithina length t the differential equation of the 
elastic line 

E J r y r " = -P (e + y r ) , 
is applicable, which through v r = e + y r reduces to: 

E J r v r " = -F v r (1) 
The general solution of this equation is as known; 

Pv r = A r coscp r +Brsincp r (2) 
with A r and B r as integration constants, while 



co r = J- ( 3 ) and k r = (4) 



For each of the n lengths an equation corresponding to 
equation (2) can be established, with 2 n temporarily 
unknown integration constants, for whose determination 
the following conditions are available: 



1) 


1 


t i me 


yi = 


0 




for 


ii = 0 






2) n - 


1 


t i me 


y r -i 




y r 


for 


x r-i = ^r-i 


and 


x r = 0 


3) n - 
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t ime 


y'r-i 




y r 


for 
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x r = 0 


4) 


1 


time 


y n = 
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or altogether 2 n -. conditions from which the values A 
and 33 may be obtained. The differentiation of equation 
(2) gives: 

k r P v r > = ~A r sin cp r + B r cos cp r (5) 

and the differentiation of equation (5) with considera- 
tion of equation (4) gives: 

k r 2 P v r » = -A r cos cp r - B r sin cp r (6) 

i.e., E J r v r " = -P v r . 

Then the determinative equations read with 

l r 

ttr = k7 (7) 

as follows; 

1. A 1 = P e 

2. A r _j cos a r _ x + B r -i sin a r _ x - A r = 0 

3. -A r _! sin a r -i + B r-i cos a r ~i - ~~~ — B r = 0 

^r 

4. A n cos an + B n sin an = P e 

The equation scheme is as follows; 

The constants A r and B r are obtained as quotients 

of two determinants and may be written in equation (2) . 
The bar then fails under a buckling load if the right-hand 
side of equation (2) becomes great beyond all limits. 
Since cos cp r and sin cp r are consistently 1, A r * and 

B r must become very great in order that P v r goes beyond 
all limits. 



Now 



Equation System 
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A r = | (8) 

Since Z as the sum of the products of k r , sin art 
cos a y , P» and e cannot go beyond all limits for finite 
values of F, e, J r , and l r , there remains as criter- 
ion for the buckling of the bar 

N >0 (9) 

No buckling occurs for values N > 0. The bar buckles at 
N = 0. Here N is the determinant of the previously- 
cited equation system and is of the (2 n) th degree for a 
bar of n different lengths. 

Examples: For tho cases n = 2 to n = 5 tho solu- 
tions of the doterminants afford the following buckling 
conditions: 

1. n = 2 

k x cot a 2 + cot cti ka ^ 0 (10) 
which, simplified, gives with \j r » = 1 - a r cot a r 

-± v s » + -f *V £ 1 (11) 

2. n = 3 

ki k 3 

k ! cota 2 cota 3 +cota,! k 2 cot a 3 +cot a x cot a 2 k 3 S 0 (12) 

k a 

For the symmetrical case l x = l 3 and J x = J 3 the ex- 
pression reduces to 

k k 

2 cot cti cot as + cot 8 ai - ~ > 0 (13) 
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S* n = 4 

k 2 k 4 

+ k ]. cot as cot cc 3 cot - cot a x 

kj 

k 

+ cot a x k 2 cot a 3 cot a.4 - cot a 



ki k 4 



k 3 

k l k 4 

+ cot a. cot a P k cot a4 - cot a, — - 

1 d 3 3 V 



k x k 3 

+ cot a x cot a a cot cc 3 k 4 - cot a 4 > 0 (14) 

kg 



4. n = 5 

+ k 1 cota 3 cota 3 cota 4 cota 5 - cota^cototg- 



+ cotai k 2 c o t a 3 co ta 4 co totg- co ta x cota. 



^2 k 5 



kj 

k 

+ cotajcotdg k 3 cota 4 cota s - cota! cota^ 



+ cotccucotagcotccg k^cotas- cota x cota 5 
+ cotaj cotctgcotiijcotai k 5 - cotctj, cota 3 

- cotog cota 4 



k S ^5 

k 2 k* 



k 3 . 



k* 
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ki ki 
- cota^cotag- 



k. 



k l k 5 

- cota, cota 



4 i 
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- cotct 3 cota 5 - 



4. k l fc 3 

- cota 4 cota 5 — — — 



k, k_ k s 

kg k4 



(15) 



The mathematical solution of equations (14) and (15) 
is already quite elaborate and unusually sensitive, be- 
cause the result appears as small difference of great 
numbers. As n increases the terms become consistently 
more extensive and so become void as far as practical 
application is concerned. 



As it is impossible to express explicitly the buck- 
ling load p contained in a and k from the preceding 
expressions, a method covering the usual cases of n = S 
and n = 3 by means of a substitute inertia moment is 
given. 



III. THE SUBSTITUTE MOMENT OF INERTIA 



Case 1 (fig. 5) . - The bar consists of n = 2 lengths 
and l 3 so that l x + l a = l, with the moments of 

and J 2 (fig. 3). Force 



I. 



inertia J 



p 

■ r n=a 



is the ulti- 



mate buckling load of the bar and complies with equation 
(10) 



k x cot cc 2 + cot a x = 0. 



Then visualize a bar of the same length I but constant 
moment of inertia J m over the entire length, whereby, 
respectively, 



J i > J m > J s and 



(16) 



J m is to be such as to comply with 
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E J m TT 2 



V-a= ^— (17) 



With the simplified relation P in place of T n=s , equa- 
tion (17) gives: 



Then 



M j 3 i /j s 

and k 2 = / = - / — (is) 



(19) 



Having recourse to equation (18) transforms equation 
(10) into: 

Ji 

^- cot a 3 = - cot a x (20) 



Herefrom the correlated cc 2 and l 1 : l 3 can be so 
determined for a constant ratio J 1 : J 2 and a variable 
a that equation (20) is fulfilled. 



Because it is 



I a /J 



and results in 



1 + € I 1 + € 
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after which: 

M m a, t 



(21) 

J 1 " t. 



.T _ ~ « 1 



(22) 

TT l g 



/ *J l 

parameter p = > — against the abscissa I,: I and the 

V Jo 



Then the values B = J can "be plotted for the 

/-— 2 



value 8 interpolated for the individual case with def- 
inite p and l x : I. Assuming J x > J 2 , S remains > 1 
and the buckling load follows at: 



E Jo TT s 



p = p 3 (23) 



The boundary curve of these curve systems follows at 
J i 

p = ~ >■ oo , wherein J 2 is finite and Ji exceeds 

all boundaries. Then it affords 



/B J x I, 

k i = — ->oo> a i = "v - >0 , and cti cot ai 

vP K i 

From the transformed equation (10) 

l x a 3 cot a 2 + l a a 1 cot a 1 = 0 

follows the condition for the boundary curve; 

I. 

— G 2 COt C 2 + 1 = 0 

U tan a I i 

A - - or _i = + 



l 2 cx 2 l v s t 
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It exists only for values 



- < a 3 < tt 



Case 2.- The "bar consists of n = 3 parts li, l 2 t 
and l lt so that l x + l 2 + l 1 ^ I and of the inertia 
mordents J 1 , J H , and J x symmetrical to the center (fig, 
4). The force P n - 3 is the buckling load. at failure and 
complies with equation (13): 

k 0 k 

2 cot an cot a 2 + cot a.7 — — — =0. 

Again a moan inertia moment J m is, introduced, so 
that: . 

S J m n 2 

Pn=3 = P = 1 (24) 

I 

As a result, equations (18) and (19) are applicable again 
and equation (13) becomes: 

777 /Ja 
2 cot cu = / — - tan a x - / — cot ai (25) 
«/ J a " 7 Ji 

From this it is possible to so determine for a con- 
stant ratio J j. : J 2 and a variable angle ax the corre- 
lated angle cc 2 and the ratio l 1 : l 3 , that equation 
(25) is complied wibh. 



For it is: 



and one obtains 
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and as before: 

£» = 5i ± and flk = 5a JL 

J x TT l x V J 2 TT l 2 

For the symmetrically stepped bar two cases 

a) J 1 > J 2 and b) J 2 > J A 

must be distinguished. 

a) (fig. 6) J x > J 2 

Pa = /— * 1 and p a = / ~ > 1 
v J 2 v Jg 

are chosen as parameters of the curves. The buckling load 
follows at 

p = 1 p a a ( 26 ) 

z 

For the boundary curve with parameter 

T 

where J 3 is finite and J x 00 we obtain 



E J I 

k i - / L >oo^ ai = — >0, and a x cot ttl >1 

P ki 

Equation (L3) may also be written as: 

2 cot cu + cot a, 8 . - -. = 0 



or + 2 cot a 2 a 2 + 1 - — - 0 a„ 2 = 0 



i 
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a s ^ i / 1 i \ 

1 + cot" a a = cot oc 2 - 2 cot a 2 — a 2 + ( — a a ) 



1 / 1 
+ = ( cot a 2 - — as 



2 ... 

d ^ 



+ = ± (cot a 2 - — a s ) 

sin a 2 ^3 ' 

1 1 + cos aa ^ 
cc 2 sin a 2 ^ a 

Only positive values of — and a P "being possible, the 

negative signs disappear. And the equation of the bound- 
ary curve becomes 



1 + cos q s _ 2j. or 00 T _ £2i 

a 2 sin a 2 l 2 a 2 l 2 



b) (fig, 7) J 2 > J i 

Here P D = / — ^ 1 is the chosen parameter. Then 
V j x 

— > 1 and the buckling load becomes 
Ji 

E Jj. TT 3 

p = _ p a (27 ) 

I 



If the parameter 

y Ji . 
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wherein J x is finite and J 3 > co , the buckling load 

of the whole bar becomes 



p k = T 



1 B Jj. TT 3 



4 I, 



2 



because the bar l 1 may for reasons of symmetry be con- 
sidered as being rigidly clamped at l 2 . For the bound- 
ary curve the condition is: 

J i _ J m n _ fl - 1 

The graphs ensuing therofrom afford the solution of the 
buckling loads Ty with the aid of the substitute moments 

of inertia, that is 

J m = J 2 p for case 1 

Jm = J 3 Pa" " " 2a 

Jm = J x P>f " " ' 2b 



IV. NOTE CONCERNING THE GRAPHS 



Charts 5 to 7 give the correlated changes in buck- 
ling load for different ratios l x : I and l 3 : I and 
different ratios of moments of inertia. They show that 
reinforcing or weakening a short end length of a compres- 
sion strut within range 

(0 < l 1 :l<0.1 or 0 < l 2 : I < 0.1) 

causes only a minor change in its buckling load, whereas, 
as figures 6 and 7, for example, show, a cross sectional 
change of a short central length alroady produces a sub- 
stantial change in the buckling load. They readily show 
the effect of local cross-sectional changes and make it 
thus possible, even for a given inertia moment and a 
stipulated buckling load P^, to determine length and 
dimension of the other piece by variation of p and Ij. : I. 
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These graphs (figs. 5 to 7) have the additional charac- 
teristic that the value p equals the parameter p in 
the extreme condition 1:1=1;' 

However, since it is impossible to include all prob- 
able cases in one graph - for l 1 : I >1 and p > 2.2 

the graphs give no values £ - the following modified 
graphs will be found to be more practical. 

Case 1 (fig. 8) . - 7 as function of l 1 : I for 
narameter 



q = 



J 2 



substitute moment of inertia: 



buckling load: 



E TT 2 



2 2 



Case 2a (fig. 9) ■ - 7 as function of 2 l 1 : I for 



the parameter 

4a 



J 2 



J x + J a 
2 



substitute inertia moment: 
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buckling load: 

E TT 2 fV+ J 2 



Case 21) (graph, fig. 10 Y . - 7^ as function of 
l 2 : I for parameter 



lb 



substitute inertia moment: 



buckling load: T , T 

„ 3 J l + J 3 
I TT ■ 



= ^ 



_2 2 



Graph. 9 loses its clearness when 2 h : I j>1. A 

close approximation in the range of 0.9 < 2 Ij. : I < 1 
is found in Fttppl 1 s solution (cf . "Lectures on Techni- 
cal Mechanics," vol. Ill, 7th ed., p. 384). 



There : 

Pk = 



TT 2 J x 



( aii + — 

j 2 ' 



V. VALIDITY IN TETMAJER'S RANGE 



Of course, it is understood that the determination 
of the failing load through the substitute moment of in- 
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ertia alone is sufficient only when the compression strut 
does not fail as a result of local over stressing. In all 
cases which fail to satisfy the previously cited condition 




the solution of the "buckling load P^ must be followed 
by further investigation, especially if the range of 
^min lies in bay center. The same applies to all struts 
with nonpr ogressively variable moment of inertia that do 
not fall in category 1, 2a, or 3b. Concerning the actual 
ultimate buckling load Pk,Br> can ^ e stated only for 

the time being that 

p k,T < p k,Br .< p k,E 

Therein the upper limit P k B is formed by that buckling 

load which satisfies the buckling condition N = 0 with 
preservation of E even in the plastic range. The lower 
limit pjj. j is obtained as follows: 



Determine 



^k,E = (> a ) 



Then let T denote Tetmajer'~s ultimate buckling stress 

at the same point A to which the stress a-^. -g belongs. 



Next , let : 

P -P 
^k,T _ ^k,E o 

k, E 

The actual buckling load P k Br then will lie be- 
tween the two limits, because, first, elastic buckling is 
ascribed to, the whole strut, while a buckling modulus 
T < E exists in the range f m i n . and secondly, the 

whole strut is given a buckling modulus T < E while T 
prevails only in the range of f m in« 
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As the buckling modulus T varies with the stress 
-only the buckling condition N ^ 0 can be investigated 
for a predetermined buckling load, whereby the buckling 
modulus T, which varies in the individual lengths l n ,' 

replaces E. The ultimate buckling load P^Br must be 

ascertained by trial, i.e., by determining the buckling 

condition IT ^ 0 for several different values 

p ( p k,T < p < p k,E^ » followed by graphical interpolation. 



Translation by J. Vanier, 
National Advisory Committee 
for Aoronautics. 
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Pigs. 1,2,3,4 




I Tetraa.jer J Euler 

r — — * 



Figure 1. 
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2 1 
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Figure 3. 



P J l J 2 

M=: J =^ 



'1 p 

=---1 -< 



l< 



Figure 4. 
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Fig. 5 




Substitute moment of inertia: J m = Jg 

Figure 5.- Buckling struts with non-progress ively variable moment 
of inertia; substitute moment of inertia J m . 
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Fig. 6 
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0 a = /— — as function of 2 1-^ '. I or Zgl I for the parameter 



Failing load Pjj- = 



E J 2 tt' 
72" 



B«3 



Substitute moment of inertia: J m = Jg J3 a ^ 



Figure 6.- Buckling struts with non-progressively variable moment 
of inertia; substitute moment of inertia J m . 
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Fig. 7 



9b - 
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^1 J'2 -Ti' "p" 



- wr/ 



■f>Pb = 2.2 



= 2.0 



= 1.8 



= 1.6 



= 1.4 



= 1.2 



= 1.0 



J m 

— as function of I 2 :I or 2 Z x : ! for the parameter 
IJ 2 




V ji 



Failing load p. = E ^JT 
Substitute moment of inertia: J m = J x 8^2 

Figure 7.- Buckling struts with non-progressively variable moment 
of inertia; substitute moment of inertia J m . 
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Fig. 



V = <J2 = 1.414- -.limit for Z x / 7— *1 and 
1.4 




V as function of l\ '. I for parameter q = 



q-*0. 
42 



Jl + J2 



Substitute moment of inertia: J m = J l+ J 2 y 2 



E 



TT 2 J l + J 2 



Buckling load P^s 



Z 



Figure 8.- Buckling struts with non-progressively variable 
moment of inertia. 
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Fig. 9 



Y a = ~ 1.414— -limit for 2^/ Z -» 1 -and *0. 
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0 .2 .4 .6 .8 1.0 
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Jl+ J2 
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Substitute moment of inertia! J m = - x V a 2 



r l+^2 



Buckling load P k = ■< ■ ■ - § . Y a 2 



Figure 9.- Buckling struts with non-progressively variable 
moment of inertia 
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Fig. 10 



Y b = V2= 1.414....1imit for Z 2 / Z -»1 anl 



!<- Z i->|<- 1 2 ->k- 1 i-H 
k- z ^ 




Vb as function of Z 2 . Z for parameter qb = 



Jl+ J2 



Substitute moment of inertia: J rn = - 

p Ji + J? 

E TT A 

Buckling load Pfc = 5 ^b^ 

Z<* 



Figure 10.- Buckling struts with non-progressively variable 
moment of inertia. 
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